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ABSTRACT 

Previous  quantitative  results  obtained  at  M.I.T.  on  the 
bending  of  steel  by  laser  line  heating  to  produce  simple 
angular  deflections  have  been  verified  and  expanded. 
Deflections  and  residual  stresses  in  compound  shapes  formed 
by  laser  line  heating  of  steel  are  measured  and  analyzed. 
From  the  measured  deflections,  algorithms  are  developed 
allowing  the  construction  of  idealized  compound  shapes 
including  the  dish,  saddle,  cone,  and  sine.  All  algorithms 
are  developed  from  single  pass  laser  line  heating  data. 
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CHAPTER  ONE 


INTRODUCTION 


It  has  been  substantially  demonstrated  by  Masubuchi  [1], 
McCarthy  [2],  et.  al.  that  bending  of  steel  plates  by  laser 
line  heating  is  practicable  for  simple  one  dimension  bending 
and  possible  for  bending  complex  shapes.  In  some  ways  laser 
bending  is  preferable  to  mechanical  forming.  These  include: 

(1)  Minimal  material  degradation  [3],  [4); 

(2)  The  ability  to  determine  deflection  angles  with  heat  in¬ 
put  parameters; 

(3)  The  ability  to  automate  the  process  fully  with  a  closed- 
loop  feedback  control  system. 

References  II]  and  [2]  also  show  that  complex  shapes  such  as 
the  dish,  saddle,  sine,  and  cone  which  are  not  readily  formed 
as  a  unit  by  mechanical  processing,  can  easily  be  formed  by 
laser  line  heating.  Much  literature  exists  extolling  the 
advantages  of  line  heating  over  mechanical  forming  for 
complex  shapes.  However,  to  date  most  attention  has  been 
focussed  on  the  use  of  flames  for  bending  steels. 

Historically,  line  heating  by  flame  has  been  an  art 
rather  than  a  science.  Attempts  have  been  made  to  quantify 
the  physics  of  flame  bending  including  the  description  of 
plastic  and  elastic  phenomena  by  boundary  integral  equations 
and  finite  element  methods  [5].  In  spite  of  the  application 


of  advanced  and  sophisticated  mathematical  tools  to  the  field 
of  flame  and  laser  bending,  simple  and  easily  programmable 
algorithms  to  predict  the  curvature  of  complex  shapes  by 
flame  or  laser  bending  do  not  exist. 

Lasers  have  been  used  for  decades  to  produce  an  intense 
and  concentrated  line  or  point  heat  for  welding  or  producing 
micro  and  mini-holes.  The  technology  already  exists  to  con¬ 
trol  accurately  the  heat  input  and  heating  pattern  of  a  laser 
by  optically  defocussing  the  beam.  This  high  degree  of 
control  indicates  the  major  advantage  of  laser  line  heating 
over  flame  line  heating.  Specifically,  a  laser  system  is,  by 
its  very  nature,  quite  adaptable  to  a  closed-loop  control 
system  for  thermal  bending  of  steel  into  complex  shapes 
including  those  with  well-defined  inflection  points. 

Today  is  a  period  of  transition  in  manufacturing  from 
simple  open-loop  and  closed-loop  control  systems  to  the 
development  and  implementation  of  "expert"  and  robotic 
systems.  The  ability  to  control  and  to  predict  accurately 
relevant  thermal  parameters  is  critical  if  plate  bending  by 
line  heating  is  to  become  an  accepted  automated  manufacturing 
process . 

In  spite  of  the  expense  of  a  high  power  laser  system,  its 
use  in  industry  is  increasing  rapidly.  Discussions  with  one 
major  high  power  CC>2  industrial  scale  laser  manufacturer 
revealed  that,  for  that  company,  the  sale  of  large  scale 
industrial  lasers  is  increasing  about  thirty  per  cent  per 


year.  Many  companies  who  would  not  be  willing  or  able  to 
apply  the  large  capital  expenditures  required  for  a  laser 
system  might  invest  in  an  automated  flame  bending  system. 
For  that  reason,  the  algorithms  developed  in  this  paper  were 
designed  to  be  applicable  to  flame  bending  with  obvious 
adjustment  to  power-deflection  parameters. 

The  purpose  of  the  research  discussed  in  this  paper  is 
threefold.  First,  to  prove  the  reproducabi 1 ity  of  laser  line 
heating  by  predicting  and  then  reproducing  one  of  the  experi¬ 
ments  discussed  in  references  (1]  and  [2],  Secondly,  to 
measure  the  residual  stresses  found  in  one  of  the  shapes 
previously  formed  by  laser  line  heating.  Finally,  and  most 
importantly,  to  develop  simple,  linear  algorithms  which 
utilyzes  the  extensive,  simple,  one  dimensional  data  base  to 
mathematically  construct  complex  shapes.  Algorithms  for  the 
dish,  saddle,  cone,  and  sine  shape  are  separately  developed. 
They  are  empirically  formulated  from  shapes  previously  formed 
by  laser  line  heating.  These  algorithms  are  then  used  to 
predict  the  deflections  of  shapes  to  be  formed  in  subsequent 
experiments  . 
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CHAPTER  TWO 

PREVIOUS  AND  POTORE  RESEARCH 

2.1  Long-Term  Research  Goals 

A  discussion  of  future  research  is  normally  included  at 
the  end  of  a  paper.  In  this  case,  in  order  to  understand 
adequately  the  intent  and  scope  of  current  research  in  laser 
forming  of  steel  plates,  one  must  be  aware  of  the  long-term 
research  goals  of  which  the  work  reflected  in  this  paper  is 
but  one  step  in  the  process. 

The  ultimate  goal  of  current  research  in  laser  bending  of 
steel  plates  is  the  development  of  a  closed-loop  and/or 
expert  system  to  fully  automate  the  process  of  forming  simple 
and  complex  shapes  frequently  found  in  the  shipbuilding 
industry.  Although  the  focus  of  research  has  been  on  using  a 
high  power  (6-15  kW)  laser  as  the  line  heater,  it  is  antici¬ 
pated  that  the  method  will  be  readily  adaptable  to  flame 
bending.  Several  subsystems  must  be  developed  in  parallel  in 
order  to  accomplish  this: 

a)  An  extensive  data  base  for  one  dimensional  bending  must  be 
established  correlating  plate  deflections  to  heating  and 
material  parameters. 

b)  Algorithms  must  be  developed  so  that  the  system  has 
initial  values  and  heating  patterns  from  which  to  start  and 
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from  which  the  system  could  estimate  correct  curvatures  for 
complex  shapes  when  encountering  new  sets  of  requirements. 

c)  A  "real-time"  deflection  measurement  subsystem  must  be 
developed  to  give  closure  to  the  feedback  system. 

d)  Computer  programs  and  methodologies  must  be  developed  to 
control  and  interface  the  database,  algorithms,  and  the  real¬ 
time  deflection  measurements.  Also,  they  must  be  capable  of 
interfacing  with  the  bending  laser  and  movable  bedplate 
systems . 

e)  Prototype  total  system  development  and  assembly. 

References  [1]  and  12]  and  part  of  this  paper  show  that  step 
"a"  has  been  accomplished.  A  low  power  laser/video  camera 
interferometer  system  is  being  developed  at  this  time.  It 
should  accomplish  step  "c".  The  main  purpose  of  this  paper 
is  to  describe  the  effort  to  accomplish  step  "b";  namely,  the 
development  of  the  algorithms.  Subsystem  "d"  and  total 
system  integration  will  follow  from  past  and  current 
research . 

2.2  Review  of  Previous  Research 

References  [1]  and  (2]  discuss  in  detail  the  results  of 
experiments  on  the  bending  of  steel  plates  by  a  high  power 
CO2  laser.  Figure  2.2-1  is  a  sketch  of  the  basic  method  of 
laser  bending  a  plate. 
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The  steel  plate  rests  on  a  movable  bed  plate  such  that  the 
stationary  beam  irradiates  the  surface  of  the  steel  as  a  line 
pass.  Primary  deflection  is  perpendicular  to  the  laser  line 
pass.  The  resulting  deflection  angle,  ©£ ,  is  the  sum  of  ©^ 
and  ©2«  All  deflection  angles  given  in  [1]  and  [2]  are  in 
terms  of  ©f. 

Many  factors  can  affect  the  amount  of  curvature.  These 
include : 

1)  Plate  thickness  and  type  of  steel. 

2)  Heat  input  which  is,  itself,  dependent  on  the  laser  power, 
optical  path  energy  loss,  absorbance  fraction  by  the  plate  of 
the  delivered  infrared  power,  and  the  speed  of  the  plate 
underneath  the  laser  beam. 

3)  Plate  restraints  such  as  clamping  and  free  edge  effects. 

4)  Plate  cooldown  rate. 

5)  Laser  beam  spot  geometry  may  have  an  effect  by  changing 
the  power  density  delivered  to  the  plates  surface. 

Since  metals  are  highly  reflective  to  the  infrared  radiation 
(“10  microns)  emitted  from  a  CO2  laser,  the  plate  surface 
must  be  coated  with  an  absorbing  substance.  Black  spray 
paint  is  currently  being  used  for  that  purpose.  Optical 
losses  are  approximately  ten  per  cent  for  the  15  kW  laser 
used  at  the  Naval  Research  Laboratory. 

Heat  input  is  being  represented  by  the  following 
parameters : 


m 

>  V.  J 


PA/V  and  P/t-/T 

where  p  is  the  laser  power  is  kWatts  less  optical  loss,  V  is 
the  plate  speed  in  inches-per-minute  and  t  is  the  thickness 
of  the  plate  in  inches.  The  deflection  angle,  0f,  is 
correlated  to  specified  P /vAT  and/or  P/t /T.  The  deflection 
angle  has  been  measured  by  dial  gages  and  a  laser 
interferometer  ([1],I2]).  Table  2.2-1  and  figures  2.2-2  and 
2.2-3  are  extracted  from  reference  [1],  Clearly,  the  results 
represented  in  the  graphs  and  table  can  be  programmed  into  a 
computer  with  routines  for  estimating  deflections  at 
parameters  not  previously  found  by  experiment.  These  results 
with  the  additional  data  measured  and  described  later  in  this 
paper  form  the  data  base  from  which  the  algorithms  for 
complex  shapes  are  developed. 

The  free  edge  of  the  plate  where  residual  stresses  must 
vanish  have  been  shown  ([1],[2])  to  have  a  pronounced  effect 
on  the  magnitude  of  the  deflection  angle  for  a  given  set  of 
parameters.  Distortion  angle  decreases  relative  to  that 
obtained  on  the  interior  of  the  plate  at  the  edges  of  a 
plate.  However,  approximately  one  beam  diameter  (1.5  inches, 
nominal)  away  from  the  edge,  this  effect  is  small  and  will 
generally  be  neglected  in  this  paper.  For  purposes  of  the 
analyses  in  this  paper  the  effects  of  the  cooldown  rate,  edge 
restraint,  and  beam  spot  pattern  also  will  be  neglected 
unless  otherwise  specified. 
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Note : 

* 

The  amount  of  angular  distortion  after  the  plate  cooled  to 
room  temperature,  or  the  final  angular  distortion,  6^  is  shown  here. 


CHAPTER  THREE 

RESIDUAL  STRESSES  IN  COMPOUND  SHAPES 
PORMED  BY  LASER  LINE  HEATING 

3.1  Introduction 

It  is  a  well-known  phenomenon  that  any  forming  of  a 
metal  will  impose  residual  stress  in  that  metal.  The  failure 
of  a  component  is  often  abetted  if  not  caused  directly  by 
residual  stresses.  When  no  external  load  is  applied  nor  any 
external  restraint  is  present  residual  stresses  can  still  be 
present.  One  cause  of  these  residual  stresses  is  a 
differential  temperature  in  the  metal.  It  is  believed  that 
the  deflections  resulting  from  line  heating  are  caused  by  a 
combination  of  plastic  and  elastic  behavior.  Since  residual 


stresses  are  limited  by  the  yield  strength  of  the  material,  a 
measurement  of  the  residual  stresses  in  a  compound  shape 
should  indicate  the  relative  contribution  of  elastic  strain 
to  the  total  deflection.  Furthermore,  for  critical 
components  where  failure  would  cause  the  loss  of  the  system, 


minimizing  or  removing  residual  stresses  might  be  required. 
This  is  particularly  true  if  the  compound  shapes  are  to  be 
welded . 

3.2  Theoretical  Relationships 

Metals  are  susceptible  to  residual  stresses  caused  by 
uneven  plastic  strains.  Masubuchi  [6]  developed  some 
fundamental  relationships  for  analyzing  residual  stresses  in 
welds.  Assuming  perpendicular  stresses  are  not  significant 

20 


the  following  development  should  apply  to  line  heating.  It 
should  be  recognized,  however,  that  since  the  heating  from  a 
laser  is  localized  to  the  surface  of  the  plate,  some  through 
thickness  residual  stresses  probably  exist.  The  following 
derivation  is  taken  from  reference  [6]. 

a.  Total  strain  results  from  elastic  and  plastic 
components : 

£x  -  £x(e)  +  Ex(p) 

£y  =  £y(e)  +  £y(p) 
frxy  *  Vxyte)  +  ^xy(p) 

where  £x,  £y,  y^Xy  are  components  of  the  total  strain 

£x(e)»  Ey(e)'  frxy(e)  are  components  of  the  elas¬ 
tic  strain 

EX(p) i  Ey(p)»  ^xytp)  are  components  of  the  plas¬ 
tic  strain 

b.  From  Hooke's  Law 

Ex(e)  =  1/E<6X  ~  ^6y ) 

6  y (e)  =  1/E(6y  -  ^6X ) 

^cy(e)  =  ^/G^Txy^ 
where  E  is  Young's  Modulus 

G  is  the  Shear  Modulus 

V  is  Poisson's  Ratio 

c.  Stresses  must  satisfy  equilibrium  conditions: 

^  6x/9  x  +  "Sfxy^y  =  ® 

■arxy/^x  +  c!6y/c)  y  =  o 

d.  Total  strain  must  satisfy  compatibility: 


and  R ' 


Then  E'  +  R'  =  0  by  compatibility. 

If  R'  is  non-zero  (uneven  plastic  strain  distribution), 
then  elastic  (residual)  stresses  must  exist  to  satisfy 
compatibility.  Implications  of  the  above  are: 

(1)  If  plastic  strains  are  linear  with  respect  to  dis¬ 
placement,  no  residual  stresses  occur. 

(2)  Residual  stresses  cannot  be  determined  by  measuring 
the  stress  change  during  external  loading  or  unloading. 

(3)  Residual  stresses  Sx,6y»  %y  can  be  calculated  if 
€x(e)»  £y ( e ) *  ^  xy ( e )  are  determined. 

The  residual  stress  is  measured  by  measuring  the  change 
in  elastic  strain  after  cutting  a  specimen.  Residual 
stresses  can  then  be  calculated  using  the  following 
equations: 

6x  =  -E/(l-t'2)  (£x  +  ^£y) 

6y  *=  ~E/  (  1-V2  )  (£y  +  i/Ex) 

where  the  above  strains  are  experimentally  measured. 

By  measuring  the  strains,  at  the  same  location,  on  both  top 
and  bottom  of  the  specimen,  one  can  determine  the  approximate 
contribution  of  planar  and  bending  stresses  to  the  total 
residual  stress.  The  average  of  the  strains  through  the 
plate  will  give  the  planar  stresses.  The  difference  in 
strains  through  the  plate  will  give  the  contribution  from 
bending . 
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3.3  Experimental  procedure 

Twenty  strain  gages  were  attached  to  the  sine  curve 
fabricated  as  one  of  the  experiments  in  reference  f  2  J .  A 
photograph  and  parameters  used  to  make  it  is  shown  in 
Appendix  A.  Ten  strain  gages  were  placed  on  both  sides  at 
the  same  x-y  coordinate.  The  plate  is  a  12"  x  12"  x  1/2" 
mild  steel  plate.  Since  the  primary  deflection  is  perpendi¬ 
cular  to  laser  beam  path  (parallel  to  the  y  axis),  the  focus 
of  the  experiment  was  determining  the  residual  stress  in  the 
x  direction.  Two  gages  per  side  were  placed  parallel  to  the 
beam  path  to  determine  if  a  "poisson"  effect  caused  strain  in 
the  y  direction.  Figure  3.3-1  shows  the  plate  geometry  and 
gage  placement. 


Figure  3.3-1:  Placement  of  Strain  Gages  on  the  Sine  Shape 


All  gages  were  BLH  Electronics  SR-4  type  FAE- 2 5- 3 5-S6 EL. 
Each  had  a  gage  factor  of  2.02  +_  1%,  a  resistance  of  3  50  +  .5 
ohms.  K  had  a  value  of  +0.26  and  the  gage  length  was  6.35  mm 


(1/4  in)  long. 

A  Vishay  Instruments  Model  P-350A  Digital  Strain 
Indicator  with  two  SB-1  Switch  and  Balance  units  was  used  to 
measure  all  strains.  This  instrument  allows  the  measurement 
of  ten  strain  gages  per  switch  unit.  One-half,  full,  or  one- 
quarter  balanced  bridge  circuits  can  be  selected.  For  this 
experiment,  the  one-quarter  bridge  configuration  was  used. 
In  this  mode  temperature  effects  can  be  automatically  compen¬ 
sated  by  connecting  an  extra  lead  from  one  of  the  strain  gage 
connections  to  an  internal  dummy  resistor  (350  ohm)  within 
the  strain  indicator.  A  5000  micro-strain  (p  in/in)  internal 
calibrator  and  a  gage  factor  input  circuit  allows  the  compen¬ 
sation  of  the  measured  resistances  for  the  wire  resistance. 

Each  strain  gage  is  balanced  to  "zero"  strain  indication 
with  the  Switch  and  Balance  unit  prior  to  cutting  the  speci¬ 
men.  After  cutting  the  specimen,  the  main  indicator  control 
is  varied  to  rebalance  the  bridge  circuit  and  the  strain  in 
microinches-per-inch  is  read  directly  off  the  Strain  Indica¬ 
tor  unit.  Table  3.3-1  shows  the  major  specifications  of  the 
Indicator  unit  ( 7 ] . 


24 


rr.v/'T-: 

.1 


V* 


Range 


Sensitivity 


Readability 


Accuracy 


» 


_+  50,000  microinch/inch  polar¬ 
ity  reversible. 


Variable.  Null  meter  deflects 
from  zero  to  full  scale  with 
40-4000  jjin/in  at  gage  factor 
=  2. 


1 y  i  n  /  i  n  . 


_+ 0 . 1  %  of  reading  or  5  pin/in 
whichever  is  greater  for  R=120 
ohms,  GF=2.  +^0.3%  of  reading 

or  5  }iin/in,  whichever  is 
greater  for  R=120  ohms,  GF=1.5 
to  4.5. 


v, 

j 


K 


i. 


y 

y .< v*  %- 


Gage  Factor  Continuously  variable  from  0.1 

to  10.0. 

Bridge  Excitation  1.5  VRMS  at  1000  Hz  square  wave. 

Table  3.3-1:  Vishay  Instrument  Strain  Indicator  and  Balance 

and  Switch  Unit  Specifications 

3.4  Experimental  Results 

Table  3.4-1  shows  the  results  of  strain  measurements  of 
the  sine  shape.  Longitudinal  spacing  is  relative  to  the  y 
axis.  Unless  otherwise  specified  all  gages  are  transversly 
placed  relative  to  the  y  a:!.s.  Microstrain  differences  and 
averages  are  between  the  top  and  bottom  of  the  specimen. 
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DIST.  FROM 
CENTERLINE 

MICROSTRAIN 

TOP 

MICROSTRAIN 

BOTTOM 

MICROSTRAIN 

DIFFERRENCE 

MICROSTRAIN 

AVERAGE 

-4 

+  121 

+  90 

+  30 

+  106 

-3 (LONG) 

-182 

-394 

-212 

-288 

-3 (TRANS) 

-582 

-458 

-1040 

-620 

-2 

-560 

-620 

-60 

-590 

-3/8 

-245 

-468 

-223 

-357 

+  3/8 

-218 

-440 

-222 

-329 

+  2 

-136 

-288 

-152 

-212 

+3 (LONG) 

-462 

-618 

-156 

-540 

+3 (TRANS) 

-1498 

+  1359 

-2857 

-695 

+  4 

-20 

-180 

-160 

-100 

Table  3.4-1:  Experimental  Measurements  of  Micro-Strain  for 

the  Sine  Shape  Formed  by  Laser  Bending 

At  point  x  *  -3  the  ratio  of  strain  in  the  longitudinal 


to  transverse  direction  is:  -182/-582  =  0.31.  At  point  x  = 

+  3,  the  ratio  is:  -462/-1498  =  0.308.  Both  values  are 

approximately  poisson's  ratio  for  steel  (^ =  0.3)  showing 
that  bending  of  this  plate  in  the  direction  parallel  to  the 
laser  pass  is  primarily  a  "poisson"  effect.  This  concept 
will  be  extremely  important  in  the  analyses  of  the  deflec¬ 
tions  of  the  complex  shapes  and  the  subsequent  development  of 
the  algorithms  for  them.  Furthermore,  on  the  basis  of  these 
results,  the  values  of  the  elastic  strain  in  the  y  direction 
at  all  points  on  the  plate  will  be  assumed  to  be:  £y(e)  =~P 
£x(e)  during  the  residual  stress  calculations. 


The  residual  stresses  due  to  bending  and  planar  effects 


now  can  be  calculated  using  the  following  equations:. 

6X  =  -E/<1-^2)  <£x  +  2/ 6y> 

6y  =  -E/d-^2)  (£y  +  ^EX) 

Ex  (based  on  experimental  evidence) 

For  steel:  E  =  2.9  x  107  psi  and  0.3.  The  equations 

become : 

bx  =  -6.37  x  1 0 1 <E x )  where  Ex  is  the  measured 

6y  =  -1.16  x  1 0 2  (£  x )  strain  in  ytin/in. 

Table  3.4-2  shows  the  residual  stresses  calculated  from  the 
experimental  data  given  in  Table  3.4-1  using  the  above 
equations.  Figure  3.4-1  shows  the  results  graphically. 


DIST.  FROM 
CENTERLINE 

MICROSTRAIN 

DIFFERENCE 

MICROSTRAIN 

AVERAGE 

X-STRESS 
BENDING 
(KSI ) 

X-STRESS 

PLANAR 

(KSI) 

-4 

30 

106 

-1.91 

-6.76 

-3 

-212 

-288 

13 . 5 

18.4 

-2 

-60 

-590 

3.80 

37.6 

-3/8 

-223 

-357 

14.2 

22.8 

+  3/8 

-222 

-329 

14  .  2 

21.0 

+  2 

-152 

-212 

9.69 

13.5 

+  3 

-156 

-540 

9.94 

34.4 

+  4 

-160 

-100 

10.2 

6.37 

Table  3.4-2:  Bending  and  Planar  Residual  Stresses  Determined 
for  the  Sine  Shape  Formed  by  Laser  Bending 


The  bending  stress  at  x  =  +2  to  the  edge  is 

approximately  constant  since  the  edge  is  clamped  on  that 
edge.  The  bending  stress  reaches  a  local  minimum  and  the 
planar  stress  reaches  a  maximum  at  approximately  the  position 
of  maximum  magnitude  of  deflection  where  the  rate  of  change 
of  the  deflection  is  a  minimum.  Therefore,  for  this  shape, 
the  bending  residual  stress  is  a  minimum  where  the  slope  of 
the  deflection  (dz/dx)  is  a  minimum.  Likewise,  the  planar 
residual  stress  is  a  maximum  where  the  magnitude  of  the 
deflection  (|z|)  is  a  maximum.  Except  for  the  region  of 
external  clamping  of  the  plate,  the  bending  stress  is 
consistently  less  than  the  planar  stress.  The  mean  ratio 
between  bending  stress  and  planar  stress  for  the  positive 
deflection  half-cycle  (X  <,  0)  is  0.43  and  the  mean  ratio 
between  bending  stress  and  total  residual  stress  is 
approximately  0.48.  The  fact  that  the  bending  residual 
stress  represents  about  one-half  of  the  total  residual  stress 
could  indicate  that  significant  changes  in  the  magnitude  of 
the  deflection  might  occur  if  stress  relieving  is  performed 
on  the  plate.  Future  experiments  could  confirm  this 
hypothesis . 

3.5  Summary  and  Conclusions 

The  results  of  the  measurement  of  residual  stresses  in 
the  sine  shape  appear  to  indicate  the  following: 


a.  The  strain  induced  in  the  plate  in  a  direction 
parallel  to  the  laser  pass  is  approximately  poisson's  ratio 
times  the  strain  in  the  direction  perpendicular  to  the 


heating  path. 

b.  The  bending  residual  stress  is  a  minimum  where  the 
deflections  are  a  maximum  since  the  magnitude  of  the  slope  of 
the  deflection  is  a  minimum  at  that  point.  The  planar  resi¬ 
dual  stress  is  a  maximum  where  the  magnitude  of  the  deflec¬ 
tions  are  a  maximum. 

c.  The  bending  residual  stress  constitutes,  on  the 
average,  of  about  50%  of  the  total  residual  stress  but  varies 
on  the  plate  from  about  10%  to  about  80%. 

d.  Although  clamping  the  plate  appears  to  make  uniform 
the  bending  stress  distribution,  clamping  does  not  appear  to 
affect  significantly  the  planar  stress  more  than  about  one 
beam  diameter  (~1.5  inch)  away  from  the  restraint. 

e.  The  planar  stress  is  approximately  symetrical  about 
the  centerline  of  the  plate  less  the  effect  from  artificially 
restraining  one  side. 

f.  No  local  variance  was  visible  in  the  planar  stress 
due  to  multiple  laser  passes.  This  indicates  that  for  over¬ 
lapping  or  close  proximity  laser  line  passes,  the  stress  is 
distributed  about  the  plate  in  accordance  with  the  geometry 
of  the  shape. 

g.  The  highest  y  direction  residual  stress  was  observed 
to  be  about  the  yield  strength  of  the  metal  (“65  ksi). 


CHAPTER  POOR 


PHASE  II  LASER  LINE  HEATING  EXPERIMENTS 

4.1  Introduction 

This  set  of  experiments  was  part  of  the  second  phase  in 
the  sequence  whose  objective  is  the  development  of  a  closed- 
loop  controlled  process  for  laser  bending  of  steel  plates. 
The  specific  goals  of  these  experiments,  conducted  at  Naval 
Research  Laboratories,  were: 

a.  Verification  of  and  testing  the  r epr oducabi 1 i ty  of 
previously  obtained  data. 

b.  Increasing  the  data  base  by  investigating  the  behav¬ 
ior  of  3/8  and  5/8  inch  thick  mild  steel  plate. 

c.  Testing  a  simple  linear  algorithm  for  predicting  the 
deflections  for  simple  curves  as  a  function  of  plate 
thickness  and  heat  delivered  to  the  metals  surface. 

d.  Quantitatively  Investigating  the  effect  of  multiple 
laser  passes  at  different  locations. 

e.  Investigating  the  effect  of  surface  preparation  on 
the  deflections  through  multiple  pass-single  loca¬ 
tion  experiments. 

As  discussed  in  Chapter  2,  reference  (1]  discusses  the 
deflection  of  steel  plates  as  a  function  of  the  heat 
delivered  to  the  metals  surface.  The  power  delivered  to  the 
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metal's  surface  is  represented  by  the  following  parameters: 

P/-J7  or  P/ts/v 

where  P  is  the  laser  power  delivered  to  the  metal,  t  is  the 
plate  thickness,  and  V  is  the  bedplate/steel  speed  beneath 
the  laser. 

It  is  anticipated  that  even  if  the  process  of  laser 
bending  steel  plates  to  produce  complex  shapes  is  controlled 
by  a  closed-loop  system,  initial  bending  to  within  twenty  to 
thirty  per  cent  of  the  final  shape  will  be  governed  by  the 
data  base  and/or  developed  algorithms.  Final  shaping  and 
dewarping  would  be  monitored  by  the  optical  laser/video 
system.  A  major  goal  of  these  experiments  was  to  test  a 
simple,  linear  algorithm  to  predict  the  one  dimensional 
deflection  of  plates  with  thicknesses  not  previously  within 
the  data  base.  To  examine  the  repr  oducabi  1  i  ty  of  this 
process,  the  deflections  of  1/2  inch  thick  plate  were 
measured  and  compared  to  that  obtained  during  Phase  1  experi¬ 
ments  . 

Since  the  completed  system/process  may  require  multiple 
passes  at  the  same  or  different  locations  on  the  plate,  both 
multiple  pass-single  location  (MPSL)  and  single  pass-multiple 
location  (SPML)  experiments  were  conducted.  It  was  expected 
that  fundamental  question  would  be  answered  concerning  any 
interaction  of  passes  from  one  location  to  another  and  inter¬ 
actions  of  multiple  passes  at  the  same  location. 
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4.2  Experimental  Procedure 

In  order  to  minimize  the  number  of  variables  and  to 
simulate  a  probable  industrial  environment,  generated  laser 
power  was  held  constant  at  7  kW  with  a  beam  diameter  of  about 
1.5  inches.  Beam  spot  geometry  was  the  "top  hat"  discussed 
in  reference  [2],  Variation  of  the  heat  linear  density 
delivered  to  the  metal  was  accomplished  by  varying  the  plate 
speed  past  the  stationary  laser  beam.  No  complex  shapes  were 
created  in  this  series  of  experiments.  Beam  path  locations 
were  predetermined  on  the  basis  of: 

a.  Maximizing  beam  pass  separation  when  independent 
measurements  were  desired  and  controlling  the  separation 
distance  during  SPML  studies. 

b.  Obtaining  no  effects  from  the  plate's  edge  parallel 
to  the  travel  direction. 

c.  Allowable  clearances  on  the  bedplate  mechanism  with 
the  steel  plate. 

All  experiments  were  performed  with  24"  x  36"  K-TEN20CF 
steel  plates.  Each  was  spray  painted  black.  No  clamping  or 
forced  cooldown  was  done.  Deflections  were  measured  with 
three  dial  gages  of  which  two  were  placed  on  one  side  of  the 
beam  pass.  Deflections  were  taken  periodically  during 
cooling  until  no  change  was  observed  (about  20  minutes  per 
pass).  Figure  4.2-1  shows  the  arrangement  for  a  typical 


experiment . 


For  an  ideal  plate,  the  minimum  deflection  at  the  centerline 
of  the  beam  pass,  Zi,2min'  should  be  zero.  In  reality, 
deflections  parallel  to  the  beam  path  are  superimposed  on  the 
plate  resulting  in  a  uniform  (except  at  the  edges)  raising  of 
the  plate.  This  is  probably  due  to  the  longitudinal  residual 
stresses  discussed  in  Chapter  Three.  Since  only  x  direction 
deflections  are  desired,  this  overall  plate  translation  in 
the  z  direction  was  removed  by  subtracting  the  minimun  value 
which  was  found  to  be  consistently  greater  than  zero.  The 
algorithms  for  complex  shapes  developed  in  the  next  chapter 
will  include  this  effect. 

The  linear  algorithm  assumes  that  is  a  linear  fun¬ 
ction  of  P/\fv  between  any  2  adjacent,  previously  determined 
data  points.  Furthermore,  it  assumes  0f  is  a  linear  function 
of  the  thickness  for  a  given  P /  between  any  two  adjacent 
thicknesses.  Deflection  angles  for  3/8"  thick  plate  was 
predicted  using  values  of  6  f  for  1/4"  and  3/8"  plates  from 
figure  2.2-1  or  table  2.2-1.  1/2"  thick  plate  was  subjected 
to  values  of  P  /  \Jv  not  previously  within  the  data  base  al¬ 
though  they  are  shown  graphically  in  figure  2.2-1.  Values  of 
6f  were  predicted  on  the  basis  of  linearity  between 
bracketing  values  of  6f  within  the  data  base  (table  2.2-1). 
This  was  done  to  test  the  linearity  of  the  curves  within 
neighborhoods  of  the  preexisting  data  and  to  test  the  repro- 
ducability  of  the  process.  Finally,  5/8"  thick  plate  was 


'T’.iT/V  V  V  V  V  V. 


’▼'T” 


subjected  to  three  different  P/\Jv  and  the  deflection  angles 


measured.  Table  4 

.3-1  gives  the  results  of  these 

expe  r iraents 

and  the  predicted 

values 

for 

3/8" 

and  1/2"  plates.  Figures 

4.3-1  and  4.3-2  are 

graphs 

of 

these 

measurements 

super i mposed 

upon  the  results  from  Phase  I 

experiments. 

THICKNESS 

P  /sfj 

P/t\/v 

6 

e 

ACTUAL 

PRED. 

1 .28 

3.41 

0.2B 

0 . 50 

1 .49 

3.98 

0.86 

1.19 

3/8" 

1 .75 

4.67 

1 .23 

1  .  20 

2.21 

5.90 

1.66 

1 .25 

1.49 

2 . 98 

0.47 

0.42 

1/2" 

1.75 

3.50 

0.73 

0 .70 

2.21 

4.43 

1 .04 

1.02 

1.49 

2.39 

0 .38 

5/8" 

1.75 

2.80 

0.82 

2.21 

3  .  54 

0.98 

Table  4.3-1:  Test  Results  Obtained  in  Parametric  Studies 

Conducted  September  1985 
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On  the  basis  of  the  experimental  results,  one  can 
observe  : 

a.  The  close  agreement  between  the  previous  and  current 
results  for  the  1/2"  plate  indicating  the  r epr oduca bi 1 i ty  of 
the  process. 

b.  For  1/2"  plating,  a  linear  algorithm  predicting 
deflections  between  data  points  previously  obtained  is 
reasonably  accurate. 

c.  1/4"  plating  should  not  be  used  to  linearly  predict 
the  behavior  of  plates  with  a  different  thickness. 

d.  The  measured  deflection  angles  for  the  5/8"  plating 
were  consistently  higher  than  expected  with  or  without  the 
use  of  a  linear  algorithm.  For  this  reason,  it  is 
recommended  that  additional  experiments  be  performed  on  5/8" 
plating  prior  to  the  integration  of  these  results  into  the 
data  base. 

Experiments  were  performed  to  determine  the  effect  of 
multiple  passes  on  the  same  plate.  P A/v  was  held  constant  at 
1.278  and  the  plate  thickness  was  3/8".  The  results  are 
given  in  table  4.3-2.  The  cumulative  (total)  deflection 
angle  as  a  function  of  the  number  of  passes  is  shown  in 
figure  4.3-3.  The  results  show  that  for  as  little  separation 
as  one  beam  diameter  (“1.5  inches)  between  passes  the 
deflection  is  approximately  constant.  This  implies  that 
under  the  conditions  of  this  experiment,  superposition  of 
angular  deflection  is  a  reasonable  assumption  in  SPIIL  mode. 


NEAR.  PASS 
( inches ) 


T.  T 
NEAR.  END 
( inches) 


(degrees ) 


A 

“ 

10.5 

CM 

B 

20.5 

5 

.23 

C 

5.75 

10.8 

.23 

D 

5.75 

16.3 

.38 

E 

3.0 

13.5 

.33 

F 

2.88 

7.88 

.31 

G 

1.50 

7.25 

.33 

Average  Angular  Deflection  Per  Pass  is  0.3  Degrees 


Die  4.3-2:  The  Effect  of  Pass  Separation  Distance  on  the 
Angular  Deflection  of  3/8"  plate  with  P / \fv  =  1.28 


Total  Angular  Deflection  (degr 


The  last  series  of  test  in  this  sequence  indicated  that, 
without  surface  preparation  (painting)  after  each  pass  in  the 
MPSL  mode,  no  significant  deflection  occurs  after  the  initial 
pass.  It  should  be  emphasized  that  in  these  experiments  the 
plate  was  allowed  to  cool  to  ambience  after  each  pass.  Since 
the  absorption  of  “10  micron  radiation  by  a  metal  increases 
rapidly  with  increasing  plate  temperature,  multiple  passes 
at  a  single  location  might  indeed  cause  significant 
deflection  if  the  material  is  not  allowed  to  cool 
significantly  prior  to  the  next  pass.  Further  investigations 
of  the  effects  of  surface  preparation  are  recommended.  This 
should  include,  at  a  minimum,  the  use  of  alternative 
preparations  such  as  white  titanium  oxide  which  is  frequently 
used  in  laser  scattering  experiments  and  should  include  the 
MPSL  processing  without  cooling  after  each  pass. 

4.4  Summary 

The  results  of  the  series  of  Phase  II  experiments  condu¬ 
cted  at  Naval  Research  Laboratories  in  Washington,  D.C.  con¬ 
tribute  to  the  understanding  of  the  behavior  of  metals  when 
subjected  to  laser  line  heating.  They  also  provided  valuable 
data  in  the  development  of  a  closed  loop  control  system. 
Specifically : 

a.  Results  for  1/2"  plating  were  verified  showing  the 
r epr oducabi 1 ity  of  the  process. 

b.  Multiply  located  passes  on  the  same  plate  cause 
deflections  which  are  approximately  linear  with  respect  to 
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the  number  of  passes.  This  was  confirmed  using  K-TEN20CF 
with  beam  path  separation  to  as  low  as  one  beam  diameter. 

c.  If  no  surface  preparation  is  performed  and  the  plate 
is  allowed  to  cool  to  room  temperature  after  each  pass  in  the 
MPSL  mode,  no  measurable  deflection  occurs  after  the  first 
pass . 

d.  A  linear  algorithm  is  sufficient  to  predict 
deflections  for  a  given  thickness  at  power  densities  not 
within  the  data  base.  A  linear  algorithm  using  1/4"  data  for 
predicting  behavior  of  plates  of  other  thicknesses  is 
inaccurate . 

Sufficient  predictability  and  control  has  been  shown  to  give 
encouragement  that  development  of  a  closed  loop  laser  bending 
system  is  practicable  for  simple  shapes. 
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CHAPTER  FIVE 

AN  ALGORITHMIC  APPROACH  TO  PREDICTING  DEFLECTIONS 
IN  A  COMPLEX  SHAPE 

5.1  Introduction 

Very  elaborate  techniques  have  been  developed  or  applied 
to  the  prediction  of  the  behavior  of  metals  undergoing  line 
heating.  in  spite  of  this,  to  date,  there  exists  no  algo¬ 
rithm  based  on  either  engineering  physics  or  a  data  base 
which  can  synthesize  a  complex,  three  dimensional  shape. 
Recalling  the  discussion  in  Chapter  Two  about  the  overall 
goals  of  present  research  in  laser  bending,  algorithms  will 
be  essential  for  automating  a  shape  synthesizing  process 
using  line  heat.  This  is  also  probably  true  for  automating  a 
metal  straightening  (dewarping)  process.  This  is  equally 
true  for  laser  and  flame  line  heating. 

The  "purest"  algorithms  are  those  which  are  based  on  a 
thorough  understanding  of  the  physics  of  the  technique 
including  the  mechanisms  of  behavior  of  the  object  in  ques¬ 
tion.  This  approach  is  usually  possible  only  for  simple 
systems  where  there  are  few  primary  variables.  Even  then, 
the  use  of  very  powerful  computers  or  experimental  equipment 
is  often  required.  Any  algorithmic  based  or  assisted  auto¬ 
mated  technique  requiring  the  use  of  a  main  frame  computer 
will  probably  not  be  cost  effective  for  a  manufacturer  such 
as  a  shipbuilder. 
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The  other  major  approach  is  to  base  the  algorithm  on 
empiracal  data.  This  approach  is  often  preferable  to  the 
"physics"  approach  since  it  allows  the  development  of  a 
process  with  potentially  many  primary  variables  and  an  incom¬ 
plete  knowledge  of  the  science  of  the  phenomenon.  A  main 
disadvantage  of  this  approach  is  that  it  may  be  more  expen¬ 
sive  to  develop  the  algorithm  because  of  the  number  of  exper¬ 
iments  required  to  create  a  sufficiently  large  data  base,  in 
the  long  run,  for  producing  an  automated  process  for  shaping 
steel  by  thermomechanical  means,  the  later  approach  is  best 
since  it  will  reduce  the  complexity  and  cost  of  the  hardware 
in  the  complete  system.  In  fact,  it  is  not  unreasonable  to 
strive  for  a  system  whose  control  system  is  maintained  and 
driven  by  a  microcomputer.  Therefore,  the  algorithms  should 
be  as  simple  and  compact  as  possible. 

The  main  objective  of  this  development  is  not  to  produce 
the  most  elaborate,  comprehensive,  or  sophisticated  algorithm 
for  laser  bending  steel  plates  into  complex  shapes.  Rather, 
it  is  to  prove  that  simple  linearized  algorithms  can  be 
constructed  which  will  produce  with  reasonable  accuracy  those 
shapes.  Furthermore,  these  algorithms  must  be  based  on  pre¬ 
existing  information  on  simple,  single  pass  deflections  al¬ 
ready  within  the  data  base.  Although  the  data  base  for 
simple  angular  deflections  is  substantial,  the  number  of 
preexisting  complex  shapes  is  not.  Therefore,  it  is  fully 
expected  that,  as  more  experiments  are  performed  on  complex 
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shapes  and  as  more  analysis  is  performed  on  the  mechanisms  of 
bending  in  response  to  line  heating,  the  algorithms  will  be 
modified  to  become  more  general  and  more  accurate.  This 
paper  develops  one  algorithm  for  each  of  the  shapes 
considered.  These  are  the  parabaloid  (dish),  the  parabolic 
hyperboloid  (saddle),  the  cone,  and  the  antisymmetric  double 
paraboloid  ("sine"). 

5.2  General  Methodology 

The  algorithms  created  in  this  chapter  are  empirically 
based  on  only  one  sample  of  that  particular  shape.  Although 
a  strenuous  attempt  was  made  to  include  known  physical  pheno¬ 
mena,  on  occasion,  underivable  factors  had  to  be  included 
which  may  or  may  not  remain  the  same  for  the  same  shapes 
constructed  under  different  conditions.  The  same  methodology 
was  applied  to  the  different  shapes  whenever  possible. 
Linearity  was  assumed  whenever  excessive  loss  of  accuracy  did 
not  result.  Symmetry  was  determined  from  geometric  proper¬ 
ties  of  the  shapes.  A  photograph  of  each  of  the  analysed 
shapes  is  included  in  Appendix  A.  These  are  extracted  from 
reference  [2].  Also  included  with  the  photograph  is  a  sketch 
of  the  pattern  of  laser  passes  and  relevant  parameters  used 
to  construct  the  shape.  The  deflections  at  specific  points 
on  the  shape  were  measured  with  a  dial  gage.  All  data  was 
compensated  for  the  thickness  of  the  plate  and  was  normalized 
such  that  the  deflection  at  the  center  of  the  plate  (x,y  =  0) 
was  numerically  zero.  The  coordinate  system  is  cartesian 
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with  the  center  of  the  plate  the  center  of  the  coordinate 
system.  The  width  of  the  plate  in  the  "x"  direction  is 


notated  Wx  and  the  width  (height)  of  the  plate  in  the  " y " 
direction  is  notated  Wy.  Unless  otherwise  noted,  the 
"top"  of  the  plate  is  that  side  such  that  the  laser  path  is 
in  the  +y  to  -y  direction.  For  a  perfectly  continuous  pure 
geometric  shape,  the  entire  shape  can  be  generated  by  func¬ 
tions  for  the  x  =  0  and  y  =  0  contours.  Figure  5.2-1  shows 
how  the  saddle  shape  can  be  generated  with  only  the  contours 
of  x  =  0  and  y  =  0  analytically  determined.  Throughout  this 
chapter  "z(x,y)M  will  represent  the  deflection  at  point  (x,y) 
and  z(0,y)  and  z(x,0)  will  represent  the  components  of  the 
deflection  from  the  x  =  0  and  y  =  0  contours,  respectively. 


ZfO,Y)=g(Y) 


Z(X,Y)=-f=(X) ♦  g(Y)-»-  CONSTANT 


Figure  5.2-1:  3  Dimensional  Shape  Generated  from  Contour  Lines 
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The  total  deflection  at  a  point  (x,y)  is: 

z(x,y)  =  f ( x  > | y— o  +  9(y)|x  =  0  +  constant 

Since  the  curve  is  normalized  such  that  2(0,01  =  0  the 
above  constant  is  zero  and: 

z(x,y)  =  f(x)|y=0  +  g(y)|x=0  or 

z (x,y)  =  z (x,0)  +  z (0,y)  . 

Geometric  symmetry  simplifies  algorithm  generation.  For 
example,  geometric  symmetry  predicts  for  the  perfect 
bowl  and  saddle:  z(x,y)  =  z(-x,y)  *  z(x,-y)  =  z(-x,-y); 

For  the  sine  shape:  z(x,y)  =  -z(-x,y); 

For  the  cone  shape:  z(x,y)  =  z(-x,y)  only. 

The  above  relationships  are  valid  only  if  the  center  of  the 
coordinate  system  is  at  the  center  of  the  shape  as  is  assumed 
for  all  the  analysis  herein. 

The  actual  complex  shapes  were  prepared  in  such  a  way  as 
to  try  to  maintain  the  symmetry.  This  was  done  by  the  choice 
of  beam  pattern  and  the  equality  of  laser  irradiation. 
However,  when  the  shapes  were  analyzed  it  was  found  that  the 
symmetry  was  not  always  preserved.  This  could  be  due  to 
uneveness  in  surface  preparation  or  plate  thickness, 
irregularities  in  laser  power  or  bed  plate  speed,  or  unequal 
edge  effects.  These  deviations  of  the  real  plates  geometry 
from  the  ideal  geometric  shape  can  be  categorized  as: 

a.  Warping:  Lack  of  fairness  or  smoothness  for  a 

particular  z-contour  line. 

b.  Asymmetry  along  an  entire  contour  line  or  surface. 


Figure  5.3.1(a)-l,  which  will  be  discussed  in  a  later  sec¬ 
tion,  shows  examples  of  both  types  of  irregularity. 
Geometry  and  the  pattern  of  laser  passes  should  cause  us  to 
predict  that  for  a  given  y  *  k  line  the  deflection  at  +x  is 
equal  to  the  deflection  at  -x.  However,  the  deflections  in 
the  y>0,  x>_0  quadrant  appear  consistently  greater  that  the 
values  in  the  y_>0,  x_<0  quadrant.  Warpage  is  also  evident  for 
the  y  =  5  line  and  from  x  =  -3  to  x  =  -4  on  the  y  =  1  line. 
For  all  the  shapes  the  edge  effects  are  somewhat  erratic. 
The  dish  was  irradiated  in  such  a  way  to  minimize  edge 
effects.  Nonetheless,  it  still  appears  in  several  of  the 
contour  lines.  For  some  of  the  other  shapes  not  fabricated 
to  minimize  end  effects  the  edge  effects  may  or  may  not  be 
considerable . 

Due  to  these  irregularities,  the  following  general 
methodology  was  used  to  analyze  the  shapes: 

a.  The  raw  (actual)  deflections  were  normalized  such 
that  z ( 0 , 0 )  =  0 . 

b.  Using  a  fairing  process  with  splines  commonly  used 
in  naval  architecture  to  smoothe  hull  shapes,  the  contour 
lines  were  "dewarped".  This  process  was  performed  such  that 
the  deflection  was  the  same  value  at  a  point  (n,k)  approached 
from  both  a  y  *  k  and  a  x  =  n  contour. 

c.  Deflections  were  averaged  in  such  a  way  as  to  obtain 
geometric  symetry  predicted  from  the  ideal  curve. 

d.  The  dewarped,  symmetrized  shape,  called  the  model 


shape,  (Zm)/  was  then  used  to  generate  the  algorithm  (ZQ). 
Although  for  each  of  the  shapes  the  algorithmically 
determined  deflections  ( Z a )  are  compared  with  the  measured 
values,  the  goal  was  to  construct  the  modeled  (dewarped  and 
symmetrized)  deflections  (Zm)  to  within  about  10  %  accuracy. 
This  procedure  can  be  justified  by  the  fact  that  plate 
warpage  and  lack  of  complete  symmetry  always  cannot  be 
controlled  or  predicted  using  a  real  process.  Secondly, 
these  algorithms  are  designed  to  initialize  the  automated 
process.  That  is,  these  algorithms  should  be  capable  of 
obtaining  "in  the  ball  park"  deflections.  It  is  assumed  that 
the  closed  loop  feedback  system  will  be  capable  of  fine 
tuning  the  shape  which  would  include  dewarping  and 
symmetrizing  the  shape  as  required  by  the  application  of  the 
shape.  it  is  estimated  at  this  time  that  approximately  70  to 
85  %  of  the  desired  shape's  deflection  will  be  obtained  by 
use  of  the  algorithm  and/or  the  data  base.  The  remaining 
deflection  will  be  controlled  by  the  feedback  system. 

The  process  of  constructing  the  model  shape  consists  of 
using  the  single  pass  ©£  (as  a  function  of  P/"\/v)  from 
reference  [2],  the  number  and  pattern  of  laser  passes, 
generalized  effects  from  the  edges,  near-pass  coupling 
effects  and,  as  required,  correlating  factors  to  generate  the 
y  =  0  and  x  =  0  contours.  These  contours  are  then  used  to 
predict  the  deflection  at  any  given  point  within  the  plate  in 
the  manner  described  above. 
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The  factor  of  two  above  comes  from  the  definition  of  0f. 

The  distributed  deflection  slope  is  given  by: 

/\z/Aq  =  Tan(k)  where  q  is  the  generalized 
coordinate . 

For  example,  a  1/2"  thick  plate  with  width  of  24"  and 
height  of  18"  is  symmetrically  irradiated  with  20  passes  on 
one  side  to  produce  a  dish.  The  beam  power  is  6.3  kW 
including  optical  losses  and  bed  plate/steel  speed  is  12  ipm. 
Then  P/\/v  =  2.02  and  from  table  2.2-1,  0f  =  0.90.  From  the 
above  equation: 

k  =  (20)  (0.9)  (22)  (16)/(2)  (24)  (18)  =  7.33° 

Tan ( k )  =  Tan (7.33°)  »  0.1286  =  flz/flq 
The  distributed  deflection  angle  is  the  basis  of  the  algo¬ 
rithms  of  all  the  shapes.  Za(x,0)  and  za(0,y)  are  assumed  to 
be  functions  of  the  Tan(k)  parameter.  It  should  be  stated  at 
this  time  that  the  Tan(k)  parameter  is  more  accurate  than 
superposition  (n‘0^/2)  for  three  dimensional  curvature  with 
many  passes.  However,  for  simple  one  dimensional  passes  such 
as  those  described  in  chapter  four  and  references  [1]  and 
[2],  superposition  appears  to  be  more  accurate. 


5.3  Shape  Analysis  and  Algorithm  Generation 
5.3.1  Dish  Shape 

a.  Derivation  of  the  Model  and  Algorithm 

The  dish  shape  was  produced  by  irradiating  a  plate  on 
one  side  longitudinally,  transversely,  and  diagonally  to  the 
y  =  0  coordinate  line.  This  pattern  is  shown  in  Appendix  A 
with  relevant  plate  and  heating  parameters.  Table  5. 3. 1-1 
gives  the  measured  and  normalized  (raw)  deflections  and 
figures  5. 3. 1-1  to  4  are  graphs  of  the  measured  data. 

Using  the  fairing  and  symmetrizing  procedure  described 
in  the  previous  section,  model  contours  are  produced.  Table 
5. 3. 1-2  and  figure  5. 3. 1-5  show  the  model  deflections. 

Recall  that  the  distributed  deflection  slope  is  given  by: 

Az/Aq  =  tan(k)  =  tan [ n6f ( Wx- 2 ) ( Wy- 2 ) /2WxWy ] 

For  this  dish:  n  =  30;  Wx=Wy=10";  Gf  =0.7  which  gives 
k  =  6.7 2°  and  tan(k)  =  tan  6.7 2°  =  0.1178 
Since  the  geometry  of  the  dish  which  was  used  to  generate  the 
model  has  x/y  symmetry,  the  y  =  0  and  x  =  0  contours  are  the 
same  curve.  Using  this  fact  and  the  value  of  tan(k)  above, 
the  following  equations  can  be  used  to  construct  the  model 


MODEL  UNEAR  DEFLECTIONS 


Figure  5.3.1— 5:  Model  Deflections 
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Table  5. 3. 1-2*.  Model  Deflections  After  Final  Symmetrizing 


za<x,0)  =  C^/2)tan(k) | x ] 
za(0,y)  *  Of/  2 )  tan  ( k )  |  y  j 


x 


=  o.iwx 

y  =  O.lWy 


za(x,0)  =  ^/'tan(k)  |x| 

o.iwx  < 

( X |  £  0.2WX 

za (0,y)  -  if tan (k) | y | 

0.1Wy  < 

|  y |  £  0 . 2 Wy 

za(x,0)  =  tan  (k  )  [  |  x  |  -  .2Wx(l-^0] 

.2WX  <  | X |  £  . 4W 

za(0,y)  =  tan ( k ) [ | y |  -  .2Wy<l-y)] 

. 2Wy  <  | y |  £  . 4W 

za(x,0)  =  tan ( k ) [ . 7 5 | x  j  h 

h  .2W*/-.1WX] 

•4wx  <  lxl  1  *5W 

za(0,y)  =  tan (k) f . 75 jy |  ^ 

H  .  2Wj3/-  .  lWy  ] 

.4Wy  <  | y |  £  .  5W 

The  multiplicative  factor  of  .75  in  the  slope  of  the  last 
series  of  equations  accounts  for  edge  effects.  Figure  5.3.1- 
6  shows  the  algorithmic  and  model  deflections  for  |q|  =  0. 
To  find  the  deflection  at  any  point  on  the  plate  one  now  only 
has  to  add  za(x,0)  and  za(0,y)  for  the  particular  coordinate 
(x,y)  . 

b.  Comparison  of  the  Algorithm,  Model,  and  Measured 
Deflections 

Table  5. 3. 1-3  gives  the  model  and  algorithmic 
deflections  and  the  deviations  in  per  cent  between  them.  The 
average  deviation  between  zm  and  za  is  about  8.5%  which  is 
within  the  goal  of  10%.  Outside  one  beam  diameter  radius  of 
the  center  where  flattening  is  greatest,  the  average 
deviation  is  about  7.5%.  Figures  5. 3. 1-7  and  5.3. 1-8  show 
the  predicted  and  model  curves  for  the  dish  shape.  Figures 
5.3. 1-9  and  5.3.1-10  show  the  predicted  and  measured  curves. 


Figures  5.3.1-11  and  5.3.1-12  show  a  "3-D"  perspective  of  the 
dish  shape.  Figure  5.3.1-11  shows  the  measured  deflections 
while  figure  5.3.1-12  show  a  3  times  magnification  of  the 
measured  deflections  to  better  show  the  shape  and  warpage. 
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Table  5. 3. 1-3:  Model  and  Algorithm  Deflections 
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Figure  5.3.1-11:  3-D  Perspective  of  the  Dish  Shape  with 
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5.3,2  Cone  Shape 

a.  Derivation  of  the  Model  and  Algorithm 

The  pattern  of  laser  passes  and  relevant  parameters  used 
to  produce  the  cone  shape  are  given  in  Appendix  A.  Table 
5. 3. 2-1  gives  the  actual,  normalized  data.  Figures  5. 3. 2-1 
to  3  are  graphs  of  the  raw  data.  In  this  case,  imposed  in¬ 
direction  strains  are  present  due  to  both  the  poisson  effect 
and  due  to  the  non-orthogonality  of  the  beam  passes  relative 
to  the  coordinate  system.  Therefore,  longitudinal  bending 
has  a  contribution  approximately  proportional  to  the  sine  of 
the  angle  between  the  y  axis  and  the  beam  pass. 

Analysis  of  the  data  seemed  to  indicate  a  significant 
coupling  effect  between  beam  passes.  The  x-direction 
deflections  are  significantly  greater  (2.5-5.75)  than 
expected  by  either  superposition  (n0f/2)  °r  bY  the  tan(k) 
parameter.  This  significant  deviation  from  expected  behavior 
may  be  an  inertial  effect.  As  the  angle  of  the  beam  pass 
with  the  y  axis  increases,  the  closer  the  beam  is  to  the 
corner  of  the  plate.  The  much  less  inertia  of  the  plate  on 
the  outboard  side  of  the  beam  pass  should  imply  that  the 
majority  of  the  angular  deflection  from  the  beam  pass  would 
be  manifested  on  that  outboard  side.  As  the  angle  increases, 
the  larger  this  effect  would  be.  Another  possible 


contributor  to  this  "amplification"  effect  is  the  density  of 
heating  at  the  focus  of  the  beam  pattern.  Geometric  symmetry 
dictates  that,  for  the  model,  z(x,y)  =  z(-x,y). 
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Table  5. 3. 2-1:  Raw  Deflection  Data  for  the  Cone  Shape 
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Table  5. 3. 2-1:  Raw  Deflection  Data  for  the  Cone  Shape  (cont) 


Figure  5.3.2  — 3  Raw  Deflection  Data 

for  the  Cone  Shape 


No  symmetry  exists  across  the  y  =  0  line.  Table  5. 3. 2-2  and 
figure  5. 3. 2-4  and  5. 3. 2-5  give  the  symmetrized  and  dewarped 
model  deflections.  As  one  would  predict  from  the  geometry  of 
a  cone,  the  changes  in  deflection  slopes  as  a  function  of  x 
are  significantly  greater  for  the  region  of  Y  >  0  than  for 
the  region  of  Y  <  0.  For  this  reason  one  algorithm  for  each 
region  was  developed.  Since  the  plate  was  cut  as  shown  in 
Appendix  A,  the  effective  area  term  in  the  tan(k)  parameter 
was  formulated  on  the  basis  of  the  area  enclosed  by  the 
perimeter  of  the  shape  minus  2  inches  divided  by  the  area  of 
the  plate.  using  these  values  gives: 

tan ( k )  =  tan[ (7) (.7) (282)/ (2) (356.25) ]  =  .0337 
It  was  discussed  above  that  an  "amplification"  factor  was 
evident  in  the  values  of  the  deflection  slopes  of  the  model. 
One  set  of  empiracal  equations  that  satisfy  the  observed 
increases  in  deflection  slope  are: 

f  =  1  x  <  Wytan (ax )  +  BD 

f=  Wx  sin(a2)/k  Wytan(a^)+BD  <,  x  <  Wytan(a3)  +  BD 
f=  Wxsin (a3 ) /2  *  k ’V  Wytan(a3)  +  BD  i  x  <  Wx/2 

where  f  is  the  amplification  factor;  a  3 ,  a  2 ,  a  3  are  the 
angles  of  the  beam  pass  with  the  y  axis;  k  is  the  distributed 
deflection  angle;  BD  is  the  beam  diameter;  and  ^is  poisson's 


ratio.  For  the  conditions  under  which  this  cone  was  made: 


Directly  inputting  the  above  values  for  f  into  the 
equations,  the  algorithm's  equations  become: 

2a ( 0 , y )  =  0  -Wy/6  <  y  <  Wy/ 6 

2a(0,y)  =  - (1-^)  tan (k) [ |y|  -  3]  |y|  >  Wy/6 

I.  y  =  0 

za(x,0)  =  tan ( k ) x  0  <  x  £  Wx/6 

za(x,0)  =  2 . 5tan  ( k  )  [  x  -  0.1WX]  Wx/6  <  x  <.  Wx/3 

za(x,o)  =  tan  ( k  )  [  5 . 7  5x  -  4Wx/3]  Wx/3  <  x  <_  Wx/2 

II.  y  <  0 

z a ( x , 0 )  =  t an ( k ) x  0  <x  <  Wx/6 

za(x,0)  =  2 . 5tan ( k )  [ x  -  0.1WX]  Wx/6  <  x  £  Wx/3 

za(x,0)  =  tan(k)  [ (l-V)  (5.75) x-.758Wx]  Wx/3<x<Wx/2 


III.  y  >  0 

za  ( x ,  0  )  =  tan  ( k  )  x  0  <  x  <_  Wx/6 

za(x,0)  =  tan(k) [3.675X  -  2.675Wx/6]  Wx/6<x<Wx/3 

za(x,0)=tan(k) [5.75(1-  V  )  -  1.1375WV  +11. 5/^] 

(W7 

for  wx/3  <  x  <_  Wx/2 

Za  (X ,  0  )  =  tan(k)  [  5 . 75x-l .  1 3  7  5WX  +  11.5/"] 
for  x  =  Wx/2 
Za (x , 0)  =  Za (-x, 0) 

Figure  5. 3. 2-6  and  5. 3. 2-7  are  graphs  of  the  model  and 
algorithm  y  =  0  and  x  =  0  contours. 


Deflections 
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b.  Comparison  of  the  Algorithm,  Model,  and  Measured 
Deflections 


Table  5. 3. 2-3  shows  the  predicted  and  model  deflections 
and  the  deviation  between  the  two.  The  average  deviation  is 
5.2%.  Figures  5. 3. 2-8  and  5. 3. 2-9  compare  the  predicted  and 
model  curves.  Figures  5.3.2-10  and  5.3.2-11  compare  the 
predicted  and  measured  values.  Figure  5.3.2-12  shows  a  "3-D" 
perspective  of  the  measured  values. 
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Table  5. 3. 2-3:  Model  and  Algorithm  Deflections 
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Table  5. 3. 2-3:  Model  and  Algorithm  Deflections  (continued) 
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a.  Derivation  of  the  Model  and  Algorithm 

Appendix  A  shows  the  parameters  and  the  pattern  used  to 
form  the  "sine"  shape.  The  coordinate  system  used  in  the 
development  of  the  algorithm  is  180°  from  that  shown  in  the 
photograph  in  the  appendix.  Table  5. 3. 3-1  gives  the 
norlmalized,  raw  data.  Figures  5. 3. 3-1  to  5. 3. 3-4  show  the 
raw  data.  As  can  be  seen  from  the  graphs,  the  edge  along  the 
x  =  +6  line  was  clamped.  Geometric  symmetry  dictates  that 
z(x,y)  =  -z(x,y)  and  z(x,y)  =  z(x,-y).  Deflections  in  a 
direction  parallel  to  the  beam  passes  were  found.  Their 
origin  has  been  discussed  previously.  The  clamping 
complicates  the  modeling  and  algorithm  formation  since  it 
causes  violation  of  the  symmetry.  Since  all  the  shapes  have 
and  will  assume  unrestrained  plates  only  the  half-plate  which 
is  unclamped  will  be  used  for  symmetrizing. 

The  first  attempt  to  model  the  plate  focussed  to  fitting 
a  sine  curve  to  the  shapes.  This  proved  to  be  too 
inaccurate.  Instead,  the  shape  was  modeled  as  an 
antisymmetric  double  parabola.  Table  5. 3. 3-2  and  figure 
3. 5. 3-5  shows  the  modeled  data.  The  inflection  point  is 
approximately  at  x,y  =0.  Significant  coupling  exists  between 
the  two  sides  of  the  plate.  The  result  is  that  the  portion 
of  the  plate  containing  the  line  connecting  the  two  sides  of 
the  plate  (each  side  irradiated  on  opposite  sides)  has  an 
unexpected  rate  of  change  of  curvature.  This  area  acts  as  a 
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Table  5. 3. 3-1:  Raw  Deflection  Data  for  the  Sine  Shape 
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Table  5. 3. 3-2:  Model  Deflections  for  the  Sine  Shape 
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clamp.  Therefore,  for  this  reason,  the  tan(k)  parameter  must 
be  modified  by  poisson's  ratio  from  about  x  =  -BD  to  about 


x  *  +BD.  The  tan(k)  parameter  remains  applicable  to  most  of 
the  interior  of  the  plate. 

The  following  equations  form  the  heart  of  the  algorithm 
for  the  sine  shape: 

tan (k )  =  tan[ (4) (0.7) (100)/(2) (144) ]  =  .0172 

za(0,y)  =  - . 25tan (k ) |y |  |y|  <  Wy/6 

za (°,y>  =  -tan(k) [ |y  |  -  15/Wy]  wy/6  <  |y |  <  Wy/2 
za(0,y)  =  za(0,-y)  <  0 

For  the  "negative"  deflection  half-cycle: 

za  ( x ,  0 )  =  -(L/^)tan{k)x  0  <  x  <  Wx/6 

za(x,0)  *  -tan ( k ) [ x  -  Wx {V -1 ) /6  '1^  Wx/6<x< . 94Wx/4 

za(x,0)  =  tan  ( k )  [  x  -  .846WX]  .94Wx/4  <  x 

za(x,0)  =  tan  (k )  [x  -  Wv/2]  1/.  W„  <  x  <  W„/2 

r.'is-y  x  x  -  x 

For  all  regions: 

za (-x, 0)  =  -za (x, 0) 

b.  Comparison  of  the  Algorithm,  Model,  and  Measured 
Deflections 

Table  5. 3. 3-3  gives  the  results  of  the  algorithm  and  the 
model  with  their  deviations.  Note  that  this  algorithm  is 
somewhat  inaccurate  for  computing  the  values  at  the  edges. 
This  is  in  part  due  to  the  clamping  effect  discussed  earlier. 
The  average  deviation  more  than  one  beam  diameter  from  the 
edges  is  9.2%.  Figures  5. 3. 3-6  and  5. 3. 3-7  show  the  model 

100 


and  algorithm  deflections  for  x,y  *  0.  Figure  5. 3. 3-8  shows 
the  model  and  algorithm  deflections  for  x,y  not  equal  to 
zero.  The  +/-  y  contour  deflections  for  a  given  x  are  the 
same  value  for  both  the  algorithm  and  moJ  1.  Figures  5. 3. 3-9 
and  5.3.3-10  show  a  comparison  of  the  actual  and  algorithm 
deflections.  Since  the  x  =  +6  edge  is  clamped  which  distorts 
the  actual  deflections,  only  the  x  <  0  deflections  are 
compared.  Figures  5.3.3-11  and  5.3.3-12  are  3-D  plots  of  the 
sine  shape. 
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Table  5. 3. 3-3:  Model  &  Algorithm  Defl.  for  Sine  (cont) 
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Figure  5.3.3-11:  3-D  Perspective  of  the  Sine  Shape  with 
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5.3.4  Saddle  Shape 

a.  Derivation  of  the  Model  and  Algorithm 

Appendix  A  shows  the  parameters  and  pattern  in  which  the 
plate  was  fabricated.  Table  5. 3. 4-1  and  figures  5. 3. 4-1  to 
5. 3. 4-4  show  the  actual  normalized  deflections.  For  this 
shape,  symmetry  dictates  that  z(x,y)  =  z(-x,y)  and  z(x,y)  = 
z(x,-y).  Using  the  above  rules  of  symmetry  and  the  general 
methodology  utilyzed  for  the  other  shapes,  the  model  deflec¬ 
tion  data  was  determined  and  is  given  in  table  5. 3. 4-2  and 
shown  in  figure  5. 3. 4-5. 

Ideally,  since  the  power  delivered  to  each  side  of  the 
plate  is  the  same,  the  magnitudes  of  the  deflections  should 
be  the  same.  This  is  not  the  case.  The  magnitudes  of  the 
deflections  of  the  X  =  0  contour  are  much  more  than  predicted 
from  the  tan(k)  parameter  or  superposition.  It  is  believed 
that  the  poisson  effect  which  has  been  discussed  is 
responsible.  The  "top"  plate  irradiation  imposes  longitudi¬ 
nal  deflections  which  enhance  the  deflections  during  the 
"bottom"  plate  irradiation.  Therefore,  the  side  which  has 
been  irradiated  second  will  have  more  deflection  than 
predicted  by  superposition  or  tan(k)  for  single  side 
irradiation.  This  increase  should  be  predictable  from  pois- 
son's  ratio  and,  in  fact,  the  plate  can  be  reasonably  modeled 
over  much  of  the  area  by  assuming  a  multiplicative  (1  +  'll ) 


factor  with  tan(k).  It  appears  from  the  data  that  the  second 
side  passes  do  not  reduce  the  first  side  (y  *  0  contour) 
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0.012 

2 

-0.376 

-0.241 

-0.136 

-0.036 

0.023 

0.044 

3 

-0.315 

-0.195 

-0.082 

0.005 

0.070 

0.095 

4 

-0.240 

-0.121 

-0.010 

0.077 

0.117 

0.156 

5 

-0.186 

-0.065 

0.021 

0.104 

0.157 

0.215 

X 

Y*0 

DEFL 

ACTUAL 

Y=1 

DEFL 

ACTUAL 

Y  =  2 
DEFL 
ACTUAL 

Y*3 

DEFL 

ACTUAL 

Y  =  4 
DEFL 
ACTUAL 

Y  =  5 
DEFL 
ACTUAL 

-5 

0.237 

0.227 

0.167 

0.137 

0.023 

-0 . 068 

-4 

0.161 

0.147 

0.126 

0.037 

-0.066 

-0.145 

-  3 

0.117 

0.085 

0.055 

-0.031 

-0.139 

-0.225 

-2 

0.039 

0.009 

-0.023 

-0.090 

-0.191 

-0.303 

-1 

0.004 

-0.015 

-0.057 

-0.132 

-0.241 

-0.358 

0 

0.000 

-0.018 

-0.057 

-0.132 

-0.242 

-0.371 

1 

0.012 

-0.003 

-0.043 

-0.111 

-0.231 

-0.352 

2 

0.044 

0.036 

-0.001 

-0.071 

-0.181 

-0.298 

3 

0.095 

0.102 

0.067 

-0.002 

-0.107 

-0.235 

4 

0.156 

0.155 

0.150 

0.085 

-0.030 

-0.132 

5 

0.215 

0.214 

0.192 

0.165 

0.032 

-0.075 

Table  5. 3.4-1:  Raw  Deflection  Data  for  the  Saddle  Shape 
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DEFL 
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i 

-0.379 
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3 
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DEFL 
MODEL 
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0.165 
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3 

0.106 
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0.033 

2 

0.042 

0.025 

-0.024 

1 

0 .008 

-0.009 

-0.061 

0 

0.000 

-0.020 

-0.068 

1 

0.008 

-0.009 

-0.061 

2 

0.042 

0.025 

-0.024 

3 

0.106 

0.088 

0.033 
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0.165 

0 .147 
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Y*  -  2 
DEFL 
MODEL 

Y  *  - 1 
DEFL 
MODEL 

Y*=0 

DEFL 

MODEL 

0.186 

0.215 

0.226 

0.102 

0.147 

0.165 

0.033 

0.088 

0.106 

-0.024 

0.025 

0.042 

-0.061 

-0.009 

0.008 

-0.068 

-0.020 

0 .000 

-0.061 

-0.009 

0.008 

-0.024 

0.025 

0.042 

0.033 

0.088 

0.106 

0.102 

0.147 

0.165 

0.186 

0.215 

0.226 

Y=  3 
DEFL 
MODEL 

Y*=  4 

DEFL 

MODEL 

Y=  5 
DEFL 
MODEL 

0.132 

0 .011 

-0.123 

0.042 

-0.072 

-0.190 

-0.038 

-0.142 

-0.260 

-0.101 

-0.200 

-0.332 

-0.140 

-0.244 

-0.379 

-0.149 

-0.258 

-0.387 

-0.140 

-0.244 

-0.379 

-0.101 

-0.200 

-0.332 

-0.038 

-0.142 

-0.260 

0.042 

-0.072 

-0.190 

0.132 

0.011 

-0.123 

Table  5. 3. 4-2:  Model  Deflections  for  the  Saddle  Shape 
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deflections  which  are  reasonably  modeled  by  the  tan(k) 
parameter.  Perhaps,  the  reason  is  that  the  process  has 
insufficient  power  to  overcome  the  high  longitudinal  residual 
stresses  imposed  during  the  "top"  plate  irradiation.  On  the 
other  hand  no  such  stresses  oppose  the  longitudinal  bending 
during  the  first  series  of  passes. 

The  algorithm  equations  become: 


tan ( k )  =  tan[ (16) (0.7) (,64)/2]  =0.0626 


za(x,0)  =  .  5 [ . 33tan (k ) ] | x | 
za(0,y)  =  . 5 [ (  ) * 5tan ( k) ] | y | 

za (x , 0 )  =  .  33tan (k ) | x | 
za(0,y)  =  Cl/) *5tan(k) [y| 


! x |  =  o.iwx 
I y I  -  o.iwy 

0.1  <  | x |  <  0.2WX 
0.1  <  Jy |  <  o . 2Wy 


Za(x,0)  =  tan ( k ) [ | x |  -  .4Wx/3]  . 2WX  <  |x|  <  . 5WX 

za(0,y)  =  (1  +^)tan(k)  [  |y  |  -  .2Wy(l  -  )]  for 

• 2Wy  <  | y I  <  .4Wy 

za(0,y)  =  (V)  ,5tan(k)  l  |y  |  +  .2Wy*C]  where 
C  =  Cl/)’5  +  (if)1,5  +  -  2  for 

.4Wy  <  | y |  <  .5Wy 

za(x,y)  =  z  a ( - x , y ) ;  za(x,y)  =  za(x,-y) 


A  % 


b.  Comparison  of  the  Algorithm,  Model,  and  Measured 
Deflections 


Table  5. 3. 4-3  shows  the  predicted  and  model  deflections. 
Average  deviation  is  10.6%.  The  average  deviation  on  the 
plate  greater  than  a  distance  of  about  2/3  BD  from  the  edges 
is  9%.  Figures  5. 3. 4-6  to  5. 3.4-8  compare  the  results  of  the 
algorithm  and  the  model.  Figures  5. 3. 4-9  and  5.3.4-10 
compare  the  results  of  the  algorithm  with  the  actual 
deflections.  Figures  5.3.4-11  and  5.3.4-12  are  3-D 
perspectives  of  the  saddle  shape. 
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CHAPTER  SIX 


CONCLUSION  AND  RECOMMENDATIONS  POR  PUTURE  RESEARCH 

The  main  purposes  of  the  research  and  analyses  in  this 
paper  have  been  proven.  Namely: 

a.  The  technique  of  irradiating  a  plate  with  a  laser 
under  controlled  conditions  is  a  reproducable  and  predictable 
phenomenon.  Complex  shapes  can  be  mathematically  constructed 
with  reasonable  accuracy  using  simple,  linearized  equations 
and  well-known  aspects  of  geometric  symmetry. 

b.  Generality  between  algorithms  has  been  established 
through  the  use  of  the  TAN(k)  parameter,  poisson's  ratio, 
plate  dimensions,  and  the  same  general  methodology. 

c.  Residual  stress  analysis  appears  to  be  consistent 
with  the  general  algorithmic  development  with  respect  to 
poisson's  effect. 

d.  All  methods  used  to  establish  the  construction  algo¬ 
rithms  are  applicable  to  flame  line  heating  as  well  as  laser 
line  heating  with  appropriate  modifications  for  differences 
in  heating  parameters.  This  should  provide  sufficient  flexi¬ 
bility  for  further  development  of  a  general  line  heating 
automated  process  using  either  lasers  or  flames. 

It  was  stated  previously  that  the  algorithms  developed 
in  this  paper  are  only  a  "first  step".  Admittedly,  they  lose 
accuracy  near  the  edges  and  sometimes  near  the  points  of 
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maximum  and  minimum  curvature.  Future  research  should,  in 
part,  focus  toward  "fine-tuning"  them.  Futhermore,  the  even¬ 
tual  goal  of  the  algorithmic  development  should  be  to  genera¬ 
lize  the  shapes  into  as  few  algorithms  as  is  practicable. 
Possibly  all  conical  and  spherical  shapes  can  be  generalized 
into  a  single  algorithm.  Due  to  the  lack  of  experimental 
data,  it  was  not  possible  to  generate  general  algorithms 
which  would  cover  a  variety  of  pure  or  combined  shapes.  For 
example,  no  algorithm  was  generated  which  would  predict  with 
correlatable  accuracy  a  combined  shape  such  as  a  sine  curve 
superimposed  on  a  conical  shape.  Also,  since  no  experimental 
data  exists  on  shapes  which  were  u nsy me t r i cal ly  heated  in 
terms  of  the  P/  V  for  the  various  passes,  the  applicability 
of  these  algorithms  to  this  condition  is  untried  and  may  not 
be  accurate.  However,  futher  experimentation  including  that 
which  is  on-going  at  the  time  of  this  writing  should  clarify 
and  better  define  the  scope  of  applicability  of  these 
routines  . 

Most  research,  up  to  this  time,  has  concentrated  in 
producing  a  process  to  laser  bend  metals  without  much  atten¬ 
tion  on  the  mechanisms  and  consequences  of  the  process.  The 
practicability  of  the  process  has  been  demonstrated.  At  some 
point  the  physics  of  laser  line  heating  should  be  thoroughly 
investigated.  This  is  of  more  than  academic  importance.  For 
example,  current  studies  indicate  that  laser  beam  spot  geome¬ 
try  has  a  significant  effect  on  the  magnitude  of  deflections. 


If  this  turns  out  to  be  true,  P/  V,  in  itself,  may  no  longer 
be  the  best  parameter  to  use  to  correlate  deflection  angles. 
More  research  into  the  questions  of  power,  thermal,  and 
stress  distributions  are  recommended. 

With  the  exception  of  the  conical  shape,  all  the  complex 
shapes  formed  up  to  now  have  been  fabricated  with  beam  passes 
either  parallel  or  perpendicular  to  the  edges  of  the  plate. 
This  may  not  be  the  most  efficient  pattern.  Other  patterns 
to  produce  the  same  shape  may  reduce  the  warpage  and 
asymmetry  found  in  the  shapes  already  constructed.  A  dish 
might  be  formed  by  a  circular  or  spiralling  beam  pattern  and 
it  conceivably  could  be  more  efficient  and/or  produce  less 
warpage.  Therefore,  research  is  recommended  in  studying 
various  patterns  to  produce  a  shape  and  in  correlating  the 
total  heat  required  to  produce  the  shape  as  compared  to  the 
patterns  used  in  the  experiments  described  in  this  paper. 

Perhaps  the  most  important  of  the  recommendations  is  the 
investigation  of  using  the  laser  bending  with  closed  loop 
system  to  symetrize  and  dewarp  a  plate  made  with  the  laser. 
It  has  been  shown  that  almost  any  shape  can  be  fabricated. 
What  has  not  been  demonstrated  up  to  now  is  if  the  same  laser 
system  can  economically  produce  those  shapes  to  the  same 
tolerances  required  in  many  industries  including  ship 
building.  The  use  of  the  laser  interferometer/video  system 
will  be  invaluable  in  accomplishing  this. 
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Ongoing  research  at  the  United  Technologies  Research 
Center  (UTRC)  in  Hartford,  Connecticut  by  M.I.T.  staff  and 
students  is  based  on  increasing  the  number  of  samples  of 
complex  shapes  and,  using  the  algorithms  within  this  paper, 
to  predict  the  deflections  under  given  power  parameters.  En¬ 
closed  within  Appendix  B  are  the  predictions  for  the 
deflections  (symmetrized  and  dewarped)  of  the  shapes. 
Furthermore,  the  research  will  determine  if  "fine  tuning" 
with  a  laser  is  practicable.  That  is,  can  the  formed  plates 
be  irradiated  again  to  achieve  the  final  desired  deflections? 
All  plates  within  this  appendix  are  assumed  to  be 
18"x24"xl/2"  and  the  pattern  of  laser  passes  are  the  same  as 
for  those  plates  used  to  develop  the  algorithms.  P/vV  and 
the  number  of  laser  passes  are  varied  to  determine  the 
accuracy  of  the  algorithm. 

The  beam  spot  of  the  CC>2  laser  at  UTRC  is  annular  vice 
the  "top  hat"  spot  geometry  of  the  laser  at  NRL .  Beam 
diameter  was  varied  in  accordance  with  optical  power  theory 
to  give  approximately  the  same  power  density  as  that 
delivered  by  the  6  kW  laser  at  NRL.  The  accuracy  of  this 
approach  will  be  determined  only  after  the  analysis  of  the 
samples . 
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z  (x i 0 )  =  .  07  37 | X |  -  .2145 


8"  <  | x |  <  12" 


4"  <  |  x  |  <_  12" 
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IV.  CONE 


A.  Parameters 


-P/v/V 


2.5 


-Of 

-Number  of  Passes 


1.2° 

Same  as  Appendix 


Equations 

z(0,y)  = 

0 

M  < 

3" 

z(0,y)  = 

-0 . 0404  [ I y I  -  1.8485] 

1 Y I  > 

3" 

Y  =  ^ 

z(x,0)  = 

.0578 | x  | 

M  < 

4" 

* 

z (x,0)  = 

0.1444[  j  X |  -  2.4] 

4"  < 

1*1 

_  *  * 
<  8" 

z (x, 0)  = 

.2325] | X |  -  1.0513] 

8"  < 

M 

<  12” 

*This  equation  applies  for  all  y's  (+/-  or  0) 


*  * 


This  equation  applies  for  Y  <.  0 

Ill 

z (x , 0  )  =  .23351 | X 1  -  1.0513]  8"  <  |x|  <  12" 

Y  >  0 

z (x, 0 )  =  . 2123 [ j  X  |  -  .6181) 


4"  <  | X |  <  12" 
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C.  Results 
X 


0 

0 

0 

0 

-6 

-.12: 

3 

0 

.173 

3 

-6 

.05: 

6 

0 

.520 

6 

-6 

.39' 

9 

0 

1.14 

9 

-6 

.921 

12 

0 

2.14 

12 

-6 

1.6 

0 

3 

0 

0 

9 

-.24 

3 

3 

.173 

3 

9 

-.07 

6 

3 

.656 

6 

9 

.41 

9 

3 

1.29 

9 

9 

1.0 

12 

3 

1.93 

12 

9 

1.6 

0 

-3 

0 

0 

-9 

-.24 

3 

-3 

.173 

3 

-9 

—  .  0  7 1 

6 

-3 

.520 

6 

-9 

.27 

9 

-3 

1.04 

9 

-9 

.79 

12 

-3 

1.74 

12 

-9 

1.5 

0 

6 

-.121 

3 

6 

.052 

6 

6 

.535 

9 

6 

1.17 

12 

6 

1.81 
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V. 


Sine 

A.  Parameters 


-P/vAT 

~6f 

-Number  of  Passes 


B.  Equations 

tan (k )  =  tanl (8) (1.2) ( .815)/2] 


v.  ••  '* 
W\« 


2.5 

1.2° 

8  Passes  Each  Side 


*  0.0684 


LL_ 


z (x,0) 

— 

-0.2279X 

x  <  4" 

2 (x, 0) 

= 

-  .  068 4x  -  .6383 

4"  <  x  < 

5.4" 

2  (  X  ,  0  ) 

a 

.0684X  -  1.3885 

5.4"  <  x 

<.  7.3 

2  (  X  f  0  ) 

* 

.  19 31 4x  -  2.5108 

7.3"  <  x 

£  12" 

2  (  X  ,  0  ) 

a 

—  2  (  —X , 0 ) 

2 (0,y) 

a 

- .  05129 | y  j 

I  Y  |  4" 

?:■#.- 


> 

V/ 

•v>. .. 

■_*  -  ™ 

V 

» 


2(0, y)  *  - . 0 6 8 4  I y I  +  .0570 


4"  <  | y |  <  9" 
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A.  Parameters 


V, 


-P/Jv 

1.708 

-0f 

0.417 

-Number  of 

passes 

Same  < 

Equations 

z (x , 0 )  =  - . 0791x 

x  <  4' 

z ( x , 0 )  =  - . 0237x  - 

.7383 

4"  <  : 

z ( x , 0 )  =  . 0237x  - 

.4818 

5.4" 

z(x,0)  =  .0670x  - 

.8713 

7.3" 

z (x, 0)  =  — z ( —x ,  y ) 

z ( 0 ,y )  =  - . 0178 | y | 

|y|  < 

z (0 ,y)  =  -.0237 

4"  < 
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SADDLE 


&.  Parameters 


-p/yv 


1.708 


0.417 


-Number  of  Passes 


16  Top;  16  Bottom 


B.  Equations 


tan ( k ) 

=  tan[ (16) ( 

.417)  ( . 8 15  > / 2 ] 

SB 

z (x, 0) 

=  .008 

M 

z (0,y ) 

=  .013 

|y| 

2 (x, 0) 

=  .0158  |  x  | 

2.5 

z  ( 0  » y ) 

=  . 0260 | y | 

2" 

z (x , 0 ) 

=  .  0475  J  x 1 

-  .1520 

5" 

z (0 ,y ) 

=  .  06 17 | y  | 

-  .1286 

3.6 

z (0 ,y) 

=  .  0867 | y | 

-  .  3084 

7.2 

where : 

z  (x,0)  _>  0 

;  z  ( 0  ,  y )  <_  0  . 

x|  =  2.5" 


2.5"  <  | x |  15" 

2"  <  | y |  <  3.6" 

5"  <  |  X |  1  12" 
3.6"  <  |  y  |  <_  7 . 2  " 
7.2"  <  I v |  <  9" 


v \\vv- -\v. 


VIII.  SADDLE 

A.  Parameters 


-P/s^T 


-Number  of  Passes  Same  as  VII. 

B.  Equations 

tan ( k )  =  tan[ (16) (1.2) ( . 815 ) /2]  =  .1374 


.0229 

M 

.  0376 

|Y| 

. 0458 | x | 

2.5 

.0753 | y | 

2" 

. 1374 | x |  - 

.4397  5" 

.1786 | y |  - 

.3721  3.6 

.  2509 | y  j  - 

.8923  7.2 

z  (x,0)  >.0; 

z  (0,y)  <_  0  . 

| X |  =  2.5" 

I y I  =  2" 

2.5"  <  | x |  1 
2"  <  |  y  |  <.3. 
5"  <  |  x  |  <.12 
3.6"  <  | y |  < 
7.2"  <  I y 1  < 


-  V.  .*  /. 


